Dynamics of Activated Escape, and Its Observation in a Semiconductor Laser 
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We report a direct experimental observation and provide a theory of the distribution of trajectories along 
which a fluctuating system moves over a potential barrier in escape from a metastable state. The exper- 
imental results are obtained for a semiconductor laser with optical feedback. The distribution of paths 
displays a distinct peak, which shows how the escaping system is most likely to move. We argue that the 
specific features of this distribution may give an insight into the nature of dropout events in lasers. 
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Fluctuation-induced escape plays an important role in 
many physical phenomena, from traditionally studied dif- 
fusion in solids and protein folding to switching in lasers 
jl],D j resonantly driven trapped electrons H , and systems 
which display stochastic resonance QJ|]. In the analysis 
of escape, it is important to be able not only to calculate, 
but also to control the escape probability. To do this one 
has to know how the system moves when it escapes. 

Escape is an example of a large fluctuation. If fluctua- 
tions are small on average, for most of the time the sys- 
tem wanders near the initially occupied metastable state 
q a and only occasionally moves far away from it. The cen- 
tral idea of the theory of large fluctuations is that paths 
to a remote state qf lie within a narrow tube centered at 
an optimal path to this state q pt(t\qf, tf) where tf 
is the instant of reaching qf. Optimal paths reveal de- 
terminism of motion in large fluctuations. They can be 
observed by analyzing the prehistory probability density 
(PPD) Ph(q, t\qf,tf) for a system to have passed through 
a point q at time t provided the system had been fluctu- 
ating about the stable state for a long time and reached 
qf at time tf. For given tf — t > 0, Ph(q,t\qf,tf) should 
peak for q lying on q opt (t\qf,tf) ||. The sharply peaked 
PPDs have indeed been observed, but so far only in ana- 
log and digital simulations and for points qf lying 
inside the attraction basin of q a || . 

In the present paper we analyze the dynamics of the 
system during escape and, using a semiconductor laser 
with optical feedback, provide a direct experimental ob- 
servation of the prehistory distribution. This distribu- 
tion displays a distinct peak, as seen from Fig. 2 below. 
We show that such peak arises even for final states lying 
behind the boundary of the domain of attraction to the 
initially occupied metastable state, e.g. behind the top 
of the potential barrier in Fig. 1. We reveal qualitative 
features of the PPD, relate them to escape dynamics, and 
compare the theoretical and experimental results. 

In the analysis of escape dynamics we will use a simple 
model of an one-variable overdamped system which per- 
forms Brownian motion in a metastable potential /7(g), 



with equation of motion 

q = -u'( q ) + mm) = ^mt - o- w 

Here, £(t) is zero- mean white Gaussian noise. We as- 
sume that the noise intensity is small compared to the 
height of the potential barrier, D <C A/7, where A/7 = 
U(qb) — U(q a ), see Fig. 1 (q a and qb are the positions of the 
local minimum and maximum of U(q)). In this case the 
escape rate W oc exp(— AU/D) is small compared to the 
characteristic reciprocal relaxation time t^ 1 = U"(q a ). 
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FIG. 1. The positions of the maxima of the prehistory 
probability density with respect to the coordinate q for given 
time t (solid line), and with respect to t for given q (triangles). 
The data of simulations refer to a Brownian particle (|l]) , Q , 
the noise intensity D = 1/60 (AU/D = 10), q f = 1.2, and 
t r = 1. The dashed line shows the asymptotic results (^), ([?]) 
for t rn (q). Inset: escape from a potential well; the motion in 
the regions A-D is discussed in the text. 

On its way to a point qf behind the barrier, the sys- 
tem is expected to move differently in the four regions 
shown in Fig. 1. In the region A behind the barrier 
top it should move nearly along the noise- free trajec- 
tory q = —U' . In the region B near the barrier top, 
where \q — qb\ S Id {Id — (27)/ A) 1 / 2 is the diffusion 
length, A = \U"(qb)\), the influence of noise becomes 
substantial. The motion is diffusive and is controlled 
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by average-strength fluctuations. The system stays here 
for the Suzuki time t s = A 1 In — Qa \/l D (Hj. In the 
region C (\q — qb\ ^> Id) the system is driven by the 
noise £(t) against the regular force —U', which requires 
a strong outburst of noise. For Gaussian noise, the prob- 
abilities of different appropriate realizations of differ 
from each other exponentially strongly. Therefore there 
is an optimal realization of noise, which is much more 
probable than others. It corresponds to an optimal path 
of the system q opt ■ For fluctuations from the attractor to 
an intrawell state such path is given by |7]] 



q op t = U'(q opt ). 



(2) 



In the region D near the attractor, \q— q a \ < {Dt r ) x ^ 2 , the 
system performs small fluctuations before a large fluctu- 
ation leading to escape occurs. 

Diffusive motion near the barrier top qb gives rise to 
a strong broadening of the distribution of fluctuational 
paths. If the destination point qf approaches qb from 
inside the well, the distribution width diverges in the 
bounce-type approximation ||. As we show, the diver- 
gence disappears if one goes beyond this approximation. 
The analytic solution will be obtained assuming that 
ln[AU/D] ;§> 1. This condition is not needed for the 
physical picture of the escape dynamics to apply, as we 
demonstrate experimentally and through simulations. 

For a Markov system ([!]), the PPD can be written as 



Ph(q,t\qf,t f ) 



p{qf,tf\q h ti) 



(3) 



where p{q\, t\ \q2, £2) is the probability density of the 
transition from q2 at the instant £2 to q\ at the in- 
stant t\ (ti > £2)- We choose the initial instant ti so 
that W' 1 > tf - U > t - ti > t r . In this time 
range the system forgets its the initial intrawell state 
qi. The statistical distribution inside and outside the 
well (not too far from the barrier top) is quasistationary, 
p(q,t\qi,ti) = p(q), and can be easily calculated. 

The prehistory distribution ph has a simple form for 
q and qf lying behind the barrier top qb in the region A 
in Fig. 1. For brevity, we give Ph{q,t\qf ,tf) in the case 
where q, qf are both in the range where U(q) is parabolic 
near qb, but qf is far enough behind qb, qf — qb ^> Id- 
The transition probability density p(qf ,tf\q,t) for such 
q, qf is known 11 1, and from (0) 



p h (q,t\q f ,0) = (2zf/l D )r(q)e Xt exp [-(z - z f e Xt ) 2 ] , (4) 

(we have set tf — 0). Here, z = z(q) = [1 — 
exp(2Ai)]~ 1/2 g, Zf = z(q f ) (note that t < 0), r(q) = 
exp(q 2 ) [1 - erf(g)] /2, and q = (q - qb)/l D . 

For \t\ < A -1 , the distribution (^) has a sharp Gaussian 
peak as a function of q, with width oc Id- Behind the 
barrier, the peak lies on the noise-free trajectory q — 
U'(q) — —X(q — qb), which arrives at qj for t = 0. 



Interestingly, the PPD peak remains sharp, with width 
~ Id, even where its maximum reaches the barrier top, 
which happens for t = t(qb) = — A -1 ln[7r 1 / 2 ((7/ — qb)/ln\- 

For earlier times —t > — t(qb), the system is mostly on 
the intrawell side of the barrier, and for large \t/t{qb)\ 
the peak of ph as a function of q moves away from the 
harmonic range. Of interest is the position t m {q) of the 
peak of ph as a function of time for given q. It shows when 
the particle was most likely to pass through the point q 
before arriving at qf. Inside the well, for qb — q 3> Id, 
the time t m and the integral width of the PPD 7(g) = 
l/ph(q, t m ) are of the form 



Atr, 



ln[2(q b - q)(q f - qb)/lj}], 7 = e l<7& - q\- (5) 



From (j^), t m depends on qb — q logarithmically. In con- 
trast, 7(g) grows linearly with qb—q- It becomes paramet- 
rically larger than the distribution width 7 ~ Id oc D 1 ! 2 
at the barrier top and outside the well. 

Far from the barrier top in the region C in Fig. 1, the 
motion of the system is determined by large fluctuations 
against the force —U'(q). In this region, ph{q,t\qf,tf) 
can be obtained using the Smoluchowski equation which 
follows from (^), 

p h {q,t\q f ,tf) = J dq' p h (q,t\q',t')p h (q',t'\qf,tf). (6) 

It is convenient to choose t' in @ so that the major 
contribution to the integral over q' came from q' lying on 
the internal side of the barrier close to qb and yet away 
from the diffusion region, qb — q a ^ qb — q' 3> Id- Then 
the second integrand in (||) is given by (^). 

The distribution Ph(q, t\q' , t') as a function of q' can 
be obtained from ([?]) by solving, in the eikonal approx- 
imation, the Fokker-Planck equation for p(q',t'\q,t) — 
exp[-S{q',t' - t\q,0)/D}. To zcroth order in D, S sat- 
isfies a Hamilton-Jacobi equation for the action of an 
auxiliary dynamical system |7j]. An appropriate Hamil- 
tonian trajectory of this system gives the optimal path 
qopt{t' — t\q, 0) for a fluctuation in which the original sys- 
tem (|l|) starts at the point q and moves further away 
from the attractor fl^| . This path is given by Eq. (0). 
The major contribution to the integral (|^) comes from 
the points q' which lie close to this path. For small 
Sq' = q' — q op t(t' — t\q, 0), it suffices to keep quadratic 
in 5q' terms in S, and then Ph(q, t\q' , t') is Gaussian in 
5q' . In the appropriately chosen parameter range, the 
time t! drops out from (^|), and one obtains JT|| 



p h (q,t\q f ,0) = \X/U'(q)\Mexp(-M), 



(7) 



where M = M(q,t) = -X(q f -q b ) [q op t(-t\q,0) - q b ] /D. 

Eq. (Q) describes the distribution of trajectories along 
which the escaping system moves inside the well. This 
distribution has a distinct peak. For given q, the peak 
is located for M(q,t m (q)) = 1. From (0), the position 
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of the peak obeys the equation dt m /dq = 1/U'(q). This 
means that, inside the potential well, the particle is most 
likely to move along the optimal path (||). In a multi- 
dimensional system, the peak of ph will lie on the most 
probable escape path, which goes from the attractor to 
the saddle point . 

The distribution (Q) is strongly asymmetric, both in q 
and t. The integral width 

7 (g) = l/p h {q,t m \q f ,Q) = e\U'(q)\/X (8) 

is independent of the noise intensity and is nonmonotonic 
as a function of q. It is maximal for U"(q) — where 
the velocity along the optimal path is maximal. The 
broadening of the tube of escape paths in time comes 
largely from the area near the barrier top. However, it is 
"amplified" as it is carried away by the trajectories flow, 
and therefore it is maximal where the flow is most fast. 




FIG. 2. (color) The prehistory probability distribution of 
the scaled radiation intensity / for experimentally observed 
dropout events in a semiconductor laser. Inset: the PPD for 
a Brownian particle, obtained from simulations for the same 
parameters as in Fig. 1. 

As —t increases further, the peak of the distribution (Q) 
approaches the diffusion region \q — q a \ ~ (Dt r Y> 2 near 
the potential minimum, and the peak width (|8j) again 
shrinks down. For large —t, the PPD (|J) goes over into 
the stationary distribution p(q), which has a nearly Gaus- 
sian peak at q a with variance Dt r /2. 

We note that, in the most interesting region C, the 
positions of the maxima of ph @j with respect to q for 
given t and with respect to t for given q are different. 
This indicates that there is no well-defined most probable 
escape path in space and time, which would go from the 
metastable state all the way over the barrier top. Still one 



can tell when the escaped system passed, most probably, 
through a given point, and where the system was most 
probably located at a given time. 

The discussed qualitative features of the prehistory dis- 
tribution can be seen from the results of digital simula- 
tions for the model potential 

U(q) = q 2 /2 - g 3 /3, (9) 

The data were obtained using a standard algorithm (jl4| , 
and refer to 8000 events. 

A distinct peak of the simulated pt is seen in Fig. 2. 
The peak of ph as a function of q changes with increas- 
ing \t\ from a narrow Gaussian near qf to a broad and 
asymmetric between q a and qb , and then again to a com- 
paratively narrow Gaussian near q a . The positions of 
the peak of ph with respect to time, t m (q), and coor- 
dinate, q m (t), are compared in Fig. 1. Both curves are 
close to each other. Outside the well they practically 
coincide and closely follow the noise-free path of the sys- 
tem, dt m jdq = —1/U'(q). The motion displays a char- 
acteristic slowing down near the barrier top. Inside the 
well the peak moves close to the optimal fluctuational 
path q — U'{q). The distribution p^ becomes time- 
independent for large \t\. Therefore t m (q) is well-defined 
only for q not too close to the potential minimum q a . 

The data of simulations in Figs. I, 2 refer to the noise 
intensities D/AU = 0.1, where the asymptotic analyti- 
cal theory applies only qualitatively. In particular, the 
expressions (|j) and (|7|) for p^ in different ranges of q do 
not merge with each other smoothly, as seen from Fig. 1. 
However, there is good qualitative agreement between 
the analytical and numerical data, including the position 
of the peak and the integral width of the PPD. 

Numerical results on the standard deviation of the 
PPD a for two noise intensities are shown in Fig. 3. As 
expected, the distribution width reaches its maximum 
well inside the well, near the inflection point U"(q) = 0. 
For higher D, the maximum is less pronounced. 
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FIG. 3. Standard deviation and the position of the maxi- 
mum of the prehistory distribution at given time, for simu- 
lated Brownian motion (left panel) and a semiconductor laser 
with optical feedback (right panel) . Solid and dashed lines on 
the left panel refer to AU/D = 10 and 3, respectively. The 
scales in the panels are different (see the text). 
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The experimental observation of the prehistory distri- 
bution was made using a semiconductor laser with optical 
feedback. The setup was similar to that used before ]lf| 
and consisted of a temperature-stabilized laser diode and 
a remote flat-surface mirror. The feedback could be con- 
trolled by a variable attenuator between them. Near the 
solitary laser threshold, such a system is unstable: af- 
ter some time of nearly steady operation the radiation 
intensity drops down, then it comparatively fast recov- 
ers to the original value, then drops down again, etc. In 
the experiment, the intensity output was digitized, with 
time resolution 1 ns. To obtain the prehistory distribu- 
tion, the intensity records were superimposed backward 
in time, starting from the instant at which the intensity, 
on its way down, reached a certain level (10% above the 
extreme dropout point). The PPD obtained from 1512 
events for the feedback 15.63% is shown in Fig. 2. 

The mechanism of power dropouts is vividly discussed 
in the literature |l5|,|2|,|l6| |l8). Most authors agree that 
the role of noise in this effect is crucial. A simple model 
|l9f describes the dropouts in terms of activation escape 
of the light intensity / over a potential barrier with shape 
(^). Previous observations Jl5| ] were in agreement with 
this model, which motivated us to measure the prehistory 
distribution for dropout events. 

It is seen from Figs. 2 and 3 that the results of the ob- 
servations agree with major qualitative results on noise- 
induced escape. The experimental PPD displays a dis- 
tinct peak. The shape of this peak is similar to the shape 
of the PPD of a noise-driven system, with the light in- 
tensity / playing the role of the coordinate q (I is scaled 
by its metastable value). The peak is narrow at small 
time \t\, and displays a characteristic broadening at in- 
termediate times. For larger \t\, the peak becomes time- 
independent. From the data in Figs. 2 and 3, the relax- 
ation time of the system is t r « 2 ns. From the value 
of the escape rate W ~ 5 x 1CU 3 ns -1 found in it 
follows that, for the model (|), AU/D w 3. Using an 
estimate do ~ (D /U" '(qa)) 1 ^ 2 for a at large — t, one can 
estimate the difference 6(AJ7/£)) 1 / 2 cto in the light inten- 
sity / at the minimum and maximum of the potential 
(@). It then follows from Fig. 3b that the system goes 
through the potential maximum for t ~ —4 ns, i.e. the 
width a reaches its maximum near the potential maxi- 
mum. In combination with larger a max /ao compared to 
that in Fig. 3a for the same AU/D, this indicates that 
the model (]l|), (||) is oversimplified. However, the overall 
form of the PPD seen from the data provides an impor- 
tant argument in favor of the stochastic model of dropout 
events. We expect that it will be possible to use high- 
resolution data on the prehistory distribution in order to 
establish a quantitative model of the system. 

In conclusion, we have analyzed the dynamics of a sys- 
tem in activated escape and revealed its distinctive fea- 
tures related to the occurrence of optimal paths and to 
the motion slowing down near a barrier top. The escape 



trajectories lie within a well-defined tube, and the sys- 
tem is most likely to go through a cross-section of this 
tube at a well-defined time before it is found behind the 
barrier. For the first time, a tube of escape trajectories 
has been observed in experiment, by analyzing dropout 
events in a semiconductor laser. 
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